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Abstract 

We present a study of the topological susceptibility \ on the lattice for full QCD 
with 2 and 4 flavours of staggered fermions at zero and finite temperature T. We 
find that x presents a sharp drop across the deconfinement transition. We also study 
the dependence of x on the quark mass at T = 0: we have no conclusive evidence 
for the expected chiral behaviour. 



1 Introduction 



The topological susceptibility x plays a relevant role in understanding several low energy 
properties of QCD. It is defined as the zero-momentum two-point correlation function of 
the topological charge density Q(x) [|1], 0], 

X = Jd 4 x d„ <0|T{^(ar)Q(0)}|0) , 

Q(x) = ^ f; v {x)f« pc {x) , (i) 

where K^(x) is the Chern current 

= ^ Al(x) (d p Al{x) - ± g f abc A b p (x)A c a (x)\ , (2) 

which satisfies Q(x) = d tl K' J '(x). In these expressions g is the QCD coupling constant 
and f abc are the structure functions of SU(3). In Eq. ([[]) we keep the derivative out of the 
expectation value in the definition of x according to the continuum prescription 0. 

The topological susceptibility in the quenched theory at zero temperature has been 
widely studied. In particular, a non-zero value of x provides an explanation for the large 
mass of the rj' particle, the would-be Goldstone boson of the axial symmetry 0. In 
Ref. j|] we have found that the value of x m pure Yang-Mills theory is x — (175(5) MeV) 4 
in agreement with the phenomenological expectation . For a review of other calculations 
see || and references therein. 

In full QCD, Ward identities and current algebra relations imply in the chiral limit ^[ 

X = w (^) , (3) 

where m is the fermion mass and Nf the number of fermions. It is interesting to verify on 
the lattice the linear dependence on the mass displayed by Eq. (|3|). 

The topological susceptibility at non-zero temperature for pure Yang-Mills theory 
has also been studied. In Ref. it has been shown that it is rather constant versus T 
below the deconfinement temperature T c and it undergoes a sharp drop across the phase 
transition. It is also interesting to study the behaviour of x f° r full QCD as a function of 
T. In particular, Eq. (|3|) suggests the existence of a drop in the signal of x across T c also 
for the full QCD case |J, (see also 0). 

In the present paper we will present a lattice determination of x in full QCD for Nf = 2 
degenerate flavours of staggered fermions at finite and zero temperature and for Nf = 4 
at finite temperature. The Monte Carlo simulation was done on the APE-Quadrics in 
Milano and Pisa. 

The method used to extract x froui the lattice is described in Section 0. The results 
and technical details of the simulation for the Nf = 4 and Nf = 2 cases are shown in 
Section ^ and Section ^ respectively. The conclusions are drawn in the last section. 
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2 Method to determine % 



We have simulated the theory on a space-time lattice. The topological charge density has 
been measured by making use of the lattice operator 

1 ±4 

Q L ( X ) = Yl W^ Tr ( n ^( x ) u p*( x )) , ( 4 ) 

fJ,Up(T=±l 

where n^(x) is a plaquette in the \i-v plane and e^ vpa is a generalized Levi-Civita tensor 
which acquires an extra minus sign for each of its indices going negative. Green functions 
containing Ql{x) differ from the continuum counterparts by a finite renormalization Z ||, 
in short 

Q L (x) = Z a 4 Q{x) + 0(a 6 ) , (5) 

where a is the lattice spacing. Actually Ql{x) mixes with fermionic operators during 
renormalization but it has been shown that the off-diagonal mixings are negligible P|. 

The lattice topological susceptibility xl is defined as 

XL = ~y- , (6) 

where V is the lattice space-time volume and Ql is the total topological charge, Ql = 
Y2xQl{%)- The definition of Xl differs from the continuum expression for x i n Eq- (0) by 
contact terms which must be subtracted. The relationship between x an d Xl is 

XL = a 4 Z 2 x + M, (7) 

where M contains mixings with operators with compatible quantum numbers. 

We have measured Ql{x) after having applied two smearing steps with smearing pa- 



rameter c = 0.9 ||10fl . The smearing process allows the operator to be less sensitive to 
the quantum fluctuations. Therefore the renormalization constant Z approaches 1 and M 
approaches 0. As a consequence the statistical errors are strongly diminished ||10|| . 

We evaluate Z and M by using the "heating method" which provides a non-perturba- 



tive determination of the renormalization constants [IT]. To calculate Z, several steps of 
an updating algorithm (the same used during the Monte Carlo simulation) are applied on 
a configuration containing a charge +1 classical instanton. These updatings thermalize 
first the quantum fluctuations, responsible for the renormalization effects, and due to the 
slowing down leave the topological content unchanged. The measurement of Ql on such 
updated configurations yields Z ■ Q. As Q is known, one can extract Z. Notice that this 
procedure is equivalent to imposing the continuum value for the 1-instanton charge (in 
the MS scheme it is +1) and extracting the finite multiplicative renormalization Z by 
evaluating a matrix element of Ql. 

The additive renormalization M is obtained in a similar way. We apply a few heating 
steps with the same updating algorithm as above on a zero-field configuration. The 
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topological susceptibility is then measured. This provides the value of M, if no instantons 
have been created during the few updating steps. As explained in [|, cooling tests must 
be done to check that the background topological charge has not been changed during 
the heating. Configurations where the topological content has been changed are discarded 
from the sample. 

Once Z and M are known, we can extract a A \ from Eq. fl7|). We have fixed a in physical 
units by measuring the string tension a on a 16 4 lattice and assuming that yfa = 420 MeV. 
To this purpose, Wilson loops ranging from 1 x 1 to 8 x 8 have been evaluated by using 
smeared spatial links and a best fit has been performed for the interquark potential to the 



form V(r) — Vq + ao/r + ar [p~2 



3 Full QCD with N f = 4 



We have numerically simulated four flavours of quarks by using staggered fermions of 
bare mass am = 0.05. For the pure gauge sector of the theory the Wilson action has 
been chosen. The $-type HMC algorithm |13[ has been used for the updating. We have 
simulated at the following values of (3 = 6/g 2 =5.00, 5.02, 5.04, 5.05 5.06 and 5.10 (g is the 
lattice bare coupling). 
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Figure 1: Polyakov loop and chiral condensate at finite T for Nf = 4. 

The lattice volume was 16 3 x 4. According to Ref. |l4j], at this temporal size the 
deconfinement transition occurs at (3 C = 5.04. We have checked this value by computing 
the Polyakov loop and the chiral condensate as shown in Fig. |l[ 

In the asymmetric lattice the temperature has been determined by the definition 
T = l/(aL t ), where L t is the temporal size. 

We have checked that the distribution of topological charge is thermalized enough. In 
Fig. we show such a distribution for (5 = 5.04 after 30 cooling steps. At smaller values 
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of am (am ~ 0.01) and higher /3's, simulations are affected by a slowing down in the 



sampling 15, [L6fl. 
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Figure 2: Distribution of topological charge P(Ql) at j3 = 5.04. 

In Table [1] we list the numerical results. The lattice spacing a has been determined 
at (3 =5.00, 5.04 and 5.10 on a 16 4 lattice. For the other values of f3, it was extrapolated 
by splines. From the value at /3 C , a = 0.30(2) fm, we infer the critical temperature to be 
T c = 164(11) MeV. 



Table 1: Xi T/T c and a versus (3 for Nf = 4 on a 16 3 x 4 lattice. 



a/fm T/Tc IP' 8 x x/MeV 4 

5.00 0.31(4) 0.9677 1.61(43) 

5.02 0.306 1 0.9804 1.13(30) 

5.04 0.30(2) 1.0000 1.21(30) 

5.05 0.296 f 1.0126 0.29(15) 

5.06 0.292 1 1.0274 0.16(12) 
5.10 0.270(8) 1.1110 0.04(12) 



f Extrapolated, see text. 

The behaviour of the ratio x(T)/x(T = 0) as a function of T/T c is shown in Fig. [3]. The 
case Nf = 2 will be discussed in the next section. For Nf = 4 the value at zero temperature 
x(T = 0) is computed as the average of all results at T < T c . For comparison, we have 
included in the figure the data for the quenched theory taken from Ref . 0] . There is clear 
evidence that the topological susceptibility in full QCD with Nf = 4 presents a drop at 
T c . Apparently the drop is sharper than in the case of pure Yang-Mills. However, the 
two theories have been simulated at rather different lattice volumes (32 3 x 8 in the pure 
Yang-Mills case). Moreover, at the small values of (3 used in our present simulation, we 
expect violations of scaling in the ratio T/T c . As a consequence we prefer not to draw 
definite conclusions about the relative slopes. 
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Figure 3: X {T)/x{T = 0) versus T/T c for N f = A,N f = 2 and JV> = 0. 



4 Full QCD with N f = 2 



Full QCD with two flavours of quarks has been simulated by using staggered fermions 
of bare mass am = 0.0125 and the usual Wilson action for the pure gauge sector. The 
updating has been performed with the R-type HMC algorithm ||13|| , each trajectory con- 
sisting in 120 steps of At = 0.004 units of molecular dynamics time. We have simulated 
at /3=5.40, 5.50, 5.55, 5.60 and 5.70 both at zero and finite temperature. 



The lattice volume was 32 3 x 8 for finite T and 16 4 for T = 0. According to Ref. (T7j at 
this temporal size the deconfinement transition occurs at (3 C = 5.54(2). We have checked 
this number by studying the Polyakov loop in a similar way as in section 3. 

We have checked that the sets of configurations produced during the simulation are 
well decorrelated and the distribution of topological charge is well sampled. For the case of 
(3 = 5.70 we needed a longer simulation to obtain a well sampled distribution of topological 
charge. 

Table 2: Xi T/T c and a versus (3 for Nf = 2 on a 32 3 x 8 lattice. 





a/fm 


T/T c 


10~ 8 x x/MeV 4 


5.40 


0.1742(50) 


0.73 


5.16(1.65) 


5.50 


0.1343(72) 


0.94 


6.21(1.74) 


5.55 


0.1249(15) 


1.02 


4.75(1.69) 


5.60 


0.1076(27) 


1.18 


2.27(0.72) 


5.70 


0.0858(35) 


1.48 


1.75(1.01) 



In Table || we report the data obtained from the simulation at finite T. From the 
value extrapolated at f3 c = 5.54, we infer the critical temperature to be T c = 194(10) (15) 
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MeV where the first error is our statistical and the second one the error induced by the 
indetermination in the value of f3 c . 

In Fig. [5] the behaviour of the ratio x(T)/x(T = 0) as a function oiT/T c is shown. For 
Nf = 2 the value at zero temperature x(T = 0) has been obtained from an independent 
simulation on symmetric lattices 16 4 at the same values of (3 (see Table |3|). There is a 
clear drop in the signal when the deconfinement temperature is crossed. 

In Fig. |] we plot the topological susceptibility obtained at zero temperature for the 
five values of j3. The corresponding set of data is listed in Table |3|. Having fixed am, the 
values of the bare mass vary as shown in the upper horizontal scale. A fit with a constant 
value of the topological susceptibility gives x — (163 ± 6) 4 MeV 4 , with the statistical test 
(chi-square)/d.o.f = 0.37; a fit with a linear homogeneous dependence on m, like Eq. (g), 
gives (V''0)bare = (6.2 ± 0.8) x 10 7 MeV 3 and (chi-square)/d.o.f = 0.94. Therefore our 
present data are compatible with both behaviours, and we are not able, within our errors, 
to check Eq. (|). 

m/MeV 

14 18 20 23 29 



5.4 5.5 5.55 5.6 5.7 

P 



Figure 4: Dependence of the topological susceptibility on the bare quark mass for Nf = 2 at 
zero temperature. 



Table 3: \ versus (3 for Nf = 2 on a 16 4 lattice. 

/3 1Q- S x x/MeV 4 = 

5.40 6.25(2.11) 

5.50 8.41(2.78) 

5.55 5.96(1.99) 

5.60 9.83(3.16) 

5.70 7.19(1.67) 
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5 Conclusions 



We have studied the topological susceptibility x i n full QCD with 2 and 4 flavours of 
dynamical fermions. At zero temperature, we have not enough precision to check the 
chiral limit, Eq. (|3|). 

At finite temperature, \ stays constant for values of T up to the deconfinement tem- 
perature T c . At T c it presents a sudden drop and \ becomes compatible with zero at 
T ~ 1.5 T c for Nf = 2 and T ~ 1.2 T c for Nf = 4. This behaviour is qualitatively similar 
to that found in Ref. j| for the quenched case. 
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